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ABSTRACT 

The double-scaling limit of the supereigenvalue model is performed in the moment 
description. This description proves extremely useful for the identification of the 
multi-critical points in the space of bosonic and fermionic coupling constants. An 
iterative procedure for the calculation of higher-genus contributions to the free en- 
ergy and to the multi-loop correlators in the double-scaling limit is developed. We 
present the general structure of these quantities at genus g and give explicit results 
up to and including genus two. 
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1 Introduction 



To date the most promising discrete approach to 2D supergravity coupled to min- 
imal superconformal models is supplied by the supereigenvalue model proposed by 
Alvarez-Gaume et al. Jl]. It is formulated in terms of a collection of N Grassmann 
even and odd variables, the "supereigenvalues" , as well as a set of even and odd 
coupling constants. For a detailed review see ref. @]. There still is no progress 
towards an understanding of the model on the basis of a generalized matrix model, 
which might provide us with a geometrical picture in the form of discretized super- 
Riemann surfaces. [] Nevertheless many of the well known features of the hermi- 
tian matrix model, such as the genus expansion, the Virasoro constraints, the loop 
equations, the moment description and the loop insertion operators ||, find their 
supersymmetric counterparts in the supereigenvalue model. From this point of view 
the supereigenvalue model appears as the natural supersymmetric generalization of 
the hermitian one-matrix model. 

The supereigenvalue model is solvable nonperturbatively in coupling constants 
but perturbatively in its genus expansion. The solution is based on a set of superloop 
equations obeyed by the superloop correlators. Away from the double scaling limit 
these equations were first solved for general potentials in the planar limit in ref. 
. In the moment description the computation of higher genus contributions could 
be automatized by an iterative procedure yielding results for general potentials and 
in principle arbitrary genus ||, representing a generalization of this very effective 
method for the hermitian matrix model |4]]. Explicit results were given for genus 
one. An alternative approach was pursued by the authors of refs. [0, [8] who managed 
to directly integrate out the Grassmann-odd variables on the level of the partition 
function, uncovering the maximally quadratic dependence of the free energy on the 
fermionic coupling constants. The supereigenvalue model also displays a connection 
to supersymmetric integrable models ||. 

In order to make contact with continuum physics the supereigenvalue model has 
to be studied in its double-scaling limit. This was done for symmetric bosonic 
potentials in ref. ]l|]in the planar limit, the case of general potentials was solved 
in ref. || for genus zero, one and partially for genus two. By making use of these 
results the authors of ref. |Tl| developed a precise dictionary between continuum 



N = 1 super-Liouville amplitudes and supereigenvalue correlators. 

In this paper the methods of the iterative solution of the supereigenvalue model 
H are applied to the double-scaling limit. The moment description turns out to 
be extremely useful for the determination of the critical points in the space of cou- 
pling constants. Just as in the hermitian matrix model Kazakov multicritical points 
]TTj] appear, related to extra zeros of the eigenvalue densities accumulating at one 



endpoint of the support. We identify the scaling properties of the moments and 
basis functions, confirming the results of refs. [p], |5| . The iterative procedure may be 
optimized to only produce terms relevant in the double-scaling limit. With these 
methods at hand we are able to state the general structure of the double scaled 

1 The only attempt in this direction is ref. 0] 
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superloop correlators and the free energy. The iteration is perfectly suited for im- 
plementation on a computer algebra system. Explicit results are stated up to and 
including genus two. 

The paper is organized as follows: In section 2 we briefly review the supereigen- 
value model, the superloop equations and their iterative solution. In section 3 the 
double-scaling limit is performed. The multi-critical points in the moment de- 
scription are identified and the scaling behaviour of moments and basis functions is 
obtained. We proceed to develop the iterative process in the scaling limit, state the 
general structure of the superloop correlators and the free energy and present our 
explicit results. Section 4 finally contains the conclusions and a short discussion of 
future perspectives. 

2 Iterative Solution of the Supereigenvalue Model 

In the following we give a brief account of the iterative procedure to solve the 
supereigenvalue model genus by genus for general potentials. A more elaborate 
description may be found in ref. ||. The aim of this paper is to apply these methods 
to the double-scaling limit. 

2.1 Superloop Insertion Operators 

The supereigenvalue model [T] is defined by the partition function 

f N N 

Z = e N2p = JiYldXidOi) JliK-Xj-eiOj) e W (-Nj2[V(Xi)-0^(\i)}), (2.1) 

i=l i<j i=l 

where the integration is over a set of N bosonic and fermionic variables, denoted by 
Aj and 9i respectively. N is even. Moreover the Grassmann even and odd potentials 
take the general form 

oo oo 

V(X i ) = y £g k K k and *(Ai) = £0+1/2^, (2.2) 

k=0 k=0 

the gk and and ik+1/2 being Grassmann even and odd coupling constants, respec- 
tively. The basic observables of this model are the connected (n|m)-superloop cor- 
relators 



W{p 1 ,...,p n \q 1 ,...,q m ) = (2.3) 



n Pi - K i n Pn - K - qi - X h jm q m ~ A 



conn 

Jm 



where the expectation value is defined in the conventional way and where "conn" 
refers to the connected part. These superloop correlators act as generating func- 
tionals for higher-point correlators like (J2i ^Aj fc J2j Xj) upon expansion in and 
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qi. The connected (n|m)-superloop correlators are related to the free energy by the 
application of the superloop insertion operators 5/SV(p) and S/S^f(p): 



5 5 5 5 

W{p u ...,p n | qi,...,q m ) = x . Jlf s ••• x . Jlf s s ••• tttt — r F, (2.4) 

where 

Hence from the one-superloop correlators W( \p) and W(p\ ) (or the free en- 
ergy F) all the multi-superloop correlators can be obtained by application of the 
superloop insertion operators. 

With the normalizations chosen above, the genus expansion of the correlators 
reads 

OO 1 

W(pi,...,p n I qi,...,q m ) = Jpj W g(Pi>--->Pn I Qi, •••,9m)- (2-6) 



Similarly we have 



OO 1 

9=0 JV 



for the free energy. 



2.2 Superloop Equations 

The genus g contribution to the one-superloop correlators W( \p) and W(p\ ) may be 
found by solving the superloop equations of the model by iteration. These two equa- 
tions express the invariance of the partition function fl2.1| ) under field redefinitions 
as shown in ref. ||. 

The Grassmann-odd superloop equation takes the form 

Jc 2m p — to Jc 2m p — lo JS z 

(2.8) 

and its counterpart, the Grassmann-even superloop equation, reads 



du V'(u) i du 1 d i du $(cu) 

W( | to) + * — — —W{(j\) — -—f — — — W(w|) 



LO 



Jc 2ni p 

W(\p) 2 



Jc 2m p 

w( P | ) w'ip | ; 



LO 



2 dp Jc 2m p 
d 



LO 



p=q 



(2.9) 
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In the derivation we have assumed a one-cut structure of the loop correlators, i.e. in 
the limit iV — > oo we assume that the eigenvalues are contained in a finite interval 
[x,y]. C is a curve around the cut. Moreover eqs. (|2.8|) and ( |2.9|) encode the super- 
Virasoro constraints Gk+1/2 Z = and Lk Z = for k > — 1, which the model obeys 
by construction [[I]]. As demonstrated in ref. these equations are accessible to 
an iterative solution in genus g. The strategy consist in employing the observation 
of refs. |7|, |8| that the free energy F depends at most quadratically on fermionic 
coupling constants. By eq. ([2.4|) this allows us to write 



W(p I ) = v(p) and W(\ p) = u(p) + u(p). (2.10) 

Here v(p) is of order one, u(p) of order zero and u(p) of order two in the fermionic 
coupling constants £k+i/2- Plugging these relations and the genus expansion ( |2.6D 
into the superloop equations ( |2.8|) and ( |2.9|) yields a set of four equations at each 
genus sorted by their order in fermionic coupling constants for the quantities v g (p), 
Ug{p) and u g {p). 

The equation of order in fermionic couplings at genus g > 1 reads 

1 " ^ 1 5 

V - u (p))u g (p) = - £ uAp)u a -Ap) + 2W^p) Ug ~ 1 ^ 

-\lt^) v °- l{q) l=J (2 - n) 

where we have introduced the linear operator V defined by 

v> f(p)= I^YM m . (2 . 12) 

Jc 2m p — 00 
The equation of order 1 at g > 1 takes the form 

g-i 

(V' - u (p))v g (p) = -(§ - v (p)^jUg(p) + J2 v g'(p) u 9-g'(p) 

g'=l 

analogously the operator \1/ is defined by 

i m= ij±m M . ( 2.i4) 

JC 2,711 p — LO 

In fact for g > 1 there is no need to solve the equations of order 2 and 3, as the 
knowledge of the quantities u g (p) and v g (p) (up to a zero mode) suffices to compute 
the free energy F g and the remaining u g 0. It is the remarkable structure of eqs. 
( |2.11| ) and fl2.13| ) which allows us to develop an iterative procedure to solve for u g 



and v g . At genus g the right hand side of (|2.11|) only contains contributions from 
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lower genera, knowing these the inversion of the operator (V' — Uo(p)) then yields 
u g (p). Similarly the right hand side of eq. ( [2.13| ) contains u g and terms of lower 
genera, allowing the computation of v g (p) (up to the zero mode of the operator 
(V — uo(p)) ). Hence the seed of the solution is uq and vq. Let us now turn to an 
important tool of the solution, the change of variables from coupling constants to 
moments and the introduction of basis functions. 



2.3 Moments and Basis Functions 

The planar as well as the higher genera contributions to the loop correlators and 
the free energy are most easily expressed by introducing instead of the couplings 
{9k, ^fc+1/2} a set of Grassmann even and odd moments {M k , J k , E k , A k }. We replace 
the bosonic couplings g k by the moments [|J 



r du V'(u) 1 

= i ^-J^-^-y)]'/* - (2 - 15) 

Jk = j^P^- kTl T^w. *>l. (2.16) 

Jc 2m (oj — y) K [ (oj — x) [oj — y) \ L/Z 



and the fermionic couplings by 



duj ^{u 



c 2ixi (oj — x) 



^[(cu-x^co-y)]^, k>l 



(2.17) 



A, = / ^j^i^-xKoj-y)}^ k>l. (2.18) 
Jc 2m [oj — y) K 

One advantage of this change of variables is that the contributions at genus 
g > 1 to the loop correlators and the free energy depend only on a finite number 
of moments, e.g. F g depends at most on 2 x 3g bosonic and 2 x (3g + 1) fermionic 
moments 0. As opposed to this F g is a function of the entire set of coupling 
constants {g k , ^+1/2}- Furthermore the description in terms of moments will prove 
extremely useful in the double-scaling limit. 

For the development of the iterative procedure we further introduce the basis 
functions x^ip) an d ^^ n \p) recursively 



X 



in), 



n-1 



(P) = TT (^\P) ~ Ex (k) (p)M n . k+1 , (2.19) 



Mi 

1 n— 1 

* H (P) = ^{^(P) - J2* ik) (p)Jn-k + i), (2.20) 



fc=l 

n-1 



k=l 



where 
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€\p) = (p- *r n - 1/2 (p - y)- 1 ' 2 , ^(p) = (p - *r 1/2 (? - 2/r n ~ 1/2 



(2.21) 



X (p) and \l/( n )(p) are basis functions of the linear operator V' 



following ref. x 
appearing in the superloop equations in the sense that 



1 



(V'-u (p))x {n \p) - 
(v>-«o(p))*W(p) = (p _ y) 



(p — x) n 
1 



71 > 1, 



n > 1. 



(2.22) 
(2.23) 



Note that 0f^(p) = 4>^{p) = 4>^{p) lies in the kernel of ( V— -u (p) ) and constitutes 
the zero mode of u fl (p) mentioned above. 



2.4 The Planar Solution 

With these definitions at hand we may now state the planar (g = 0) solution of eqs. 
([O]) and O obtained in refs. §, §: 



W ( I p) 



da; y ; (a;) 

C 27TZ p — UJ 



(p-x)(p-y) 
(u - x)(u - y) 



1/2 



(2.24) 



+ ^ 52 . (5l " A ^ 4"M+ 1A2(5l " Ai: 



4 Mi(x-y) 



4 Ji{x-y) 



#(p) 



W (p 



C 27T2 p — UJ 



(u - x)(u - y) 
(p-x)(p-y) 



1/2 



Hi + Ai 



2 [(p-x)(p-j/)]V5 



(2.25) 



Here the endpoints x and y of the cut on the real axis are determined by the 
boundary conditions 



du 



c 2m J(u- x ){u-y)' 



du 
c 2m 



uV'(u) 



(u - x)(u - y) 



(2.26) 



deduced from our knowledge that W( \ p) = 1/p + 0(p~ 2 ). This is also the reason 
why there are no zero mode contributions possible for the W g { \ p). We shall make 
use of the following rewriting of the purely bosonic part of W ( | p) |12] 



Uq 



-| oo 

(p) = V\p) - - [(p - x)(p - y)} 1 ' 2 {(p - *r 1 M q + (p - yf- 1 J q }, (2.27) 

A q=l 
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derived by deforming the contour integral in eq. ( 2.24 ) into one surrounding the 
point p and the other encircling infintity. To take the residue at infinity one rewrites 
(p — uj)~ 1 as 



1 



1 



1 



p — u) 2 (p — x) — (a? — x) ' 2 (p — y) — (uo — y) ' ^ 

and expands in f^Ef) an d (i^f) respectively. Doing this for the fermionic W (p | ) 
yields 



Vo(p) = *(p) - - [(p-ar)(p- 



^ | (p _ s? + (p _ y)ff -i Ag }. (2 . 29) 



9=2 



It is important to realize that the bracketed terms in eq. ( |2.27 ) as well as in eq. 
( p.29| ) are actually identical. Here we see that the planar solution is special in the 
sense that it depends on the full set of moments. Interestingly enough this is not 
the case for higher genera. 



3 The Double— Scaling Limit 

Similar to the situation in the hermitian matrix model the "naive" iV — > oo contin- 
uum limit of the supereigenvalue model is unsatisfactory as it leaves us only with 
the planar contributions, easily seen from eqs. ( |2.6| ) and (|2.7| ). More interesting from 
the point of view of continuum physics is the double-scaling limit. Crucial for this 
scenario is the observation that there exists a subspace in the space of couplings 
{9k, ^fc+1/2} where all higher genus contributions to the free energy F g diverge. This 
enables us to take the double-scaling limit, where one simultaneously approaches 
the critical subspace of couplings and takes N — ► 00, giving contributions to the free 
energy from all genera. Let us see how this works in detail. 



3.1 Scaling of Moments and Basis Functions 

The analysis of the scaling behaviour for the bosonic quantities was carried out by 



Ambj0rn et al. P, IT2| in the framework of the hermitian matrix model. Consider the 



case of generic, i.e. non-symmetric, potentials V(p) and \P(p). The m'th multicritical 
point is reached when the eigenvalue density (uo(p) — V'(p)) of eq. ( [2.271 ), which 
under normal circumstances vanishes as a square root on both ends of its support, 
aquires (m — 1) additional zeros at one end of the cut, say x. The condition for 
being at an m'th multicritical point then simply is 



Mr = M 2 c = ---ilC_ 1 = 0, 1^0, k>m, Jf^O, 1>1, (3.1) 

defining a critical subspace in the space of bosonic couplings Denote by g\ a 
particular point in this subspace for which the eigenvalue distribution is confined to 
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the interval [x c , y c ]. If we now move away from this point the conditions of eq. ( [3.1|) 
will no longer be fulfilled and the cut will move to the intervall [x, y] . Assume we 
control this movement by the parameter a and set jp 



x = x c -aA 1/m (3.2) 
p = x c + air (3.3) 

The scaling of p and the introduction of its scaling variable 7r is necessary in order 
to speak of the double-scaling limit of the superloop correlators. A plays the role of 
the cosmo logical constant. We will now deduce the scaling of y by further assuming 
that the critical subspace {g%} is reached as 

9k = 9 ■ 9 c k , (3-4) 

where g is a function of a to be determined. Imposing the boundary conditions 
( PD yields 

y-y c ~a m , and g - 1 ~ a m . (3.5) 

Knowing this one easily computes the m'th multicritical scaling behaviour of the 
bosonic moments 

M fc ~a m - fc , fce[l,m-l], (3.6) 

while the higher M-moments and the J-moments do not scale. 

Moreover the functions (p) and 0^ (p) are found to behave like 

#(p)~a-»- x '\ ^\p)-a-'l\ (3.7) 
from which one proves the scaling behaviour of the basis functions 

X (n) (p) ~ a~ m - n+1/2 } ¥ n \p) ~ a" 1/2 , (3.8) 

following ref. ||. 

Let us now turn to the scaling of the fermionic moments 5^ and Similar 
to the bosonic case the function (vo(p) — ^(p)) of eq. ( |2.29|) usually vanishes at the 



endpoints of the cut like a square root. We will fine tune the coupling constants 
£fc + i/2 in such a manner that (n — 1) extra zeros accumulate at x, i.e. 

S c 2 = ---S^ 1 = 0, ~^0, k>n, A^O, l>2, (3.9) 

where H£ = Sfc[x c , y c , ^+1/2]- m addition the results of ref. tell us that the mo- 
ments Hi and Ai will always appear in the combination (Si — Ai). This suggests 
to impose the constraint — A^ = on these moments. As there is no analogue 
to the boundary conditions (|2.26|) for the fermionic quantities we are free to choose 
the scaling of the coupling constants £fc+i/2- We set 
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6+1/2 = a 1 ' 2 C k+1/ 2, (3-10) 

and will comment on this choice later on. From this one derives (Hi — A{) ~ a" -1 / 2 
and H/j ~ a n ~ fe+1 / 2 for A; e [2,n — 1]. All other fermionic moments scale uniformly 
with a 1 / 2 . So far the fermionic scaling is completely independent of the bosonic 
scaling, governed by the integer n. We shall, however, introduce the requirement 
that the scaling part of the bosonic one-superloop correlator of eq. ( |2.24| ) scales 
uniformly, i.e. we require (wo(j>) ~ V'ip)) an d uo(p) to scale in the same way IJ. As 
(uq(p) — V'{p)) ~ a m_1//2 we arrive at the following condition on n: 



n = m + l. (3.11) 



And therefore 



Ei - Ai ~ a m+1/2 , E fc ~ a m - k+3/2 , k e [2, m] 

(3.12) 

E fc ~ a 1/2 , k > m A t ~ a 1/2 , I > 1. 

The double-scaling limit is now defined by letting N — > oo and a — > but keeping 
the string coupling constant a = a~ 2m ~ l N~ 2 fixed. Q 

The above scaling behaviour may be shown to be equivalent to the one obtained 
in ref. []5J. Based on these results Abdalla and Zadra |TD[ proved that the central 



charge of the superconformal field theory described in the continuum is given by 

(2m -l) 2 7 

c=l-- - = 0, — m>l. (3.13) 

m 2 

Similar to the situation in the bosonic model the m — 1 fixed point may not be 
reached by the described techniques and should be treated seperately. 

3.2 The Iteration for u g (p) and v g (p) 

By making use of the above scaling properties of the moments and basis functions 
we may now develop the iterative procedure which allows us to calculate directly the 
double-scaling relevant versions of u g (p) and v g (p). The iterative scheme described 
in ref. is quite simple: By eqs. (|2.22|) and ( p.23|) every u g and v g may be written 



as a linear combination of basis functions x an d ^^ n \ where the coefficients of 
this expansion are simply read off the poles at x and y of the right hand sides of 

2 Let us now comment on the choice of eq. ( |3.10[ ). At first sight one might have expected a 
scaling like £,k+i/2 = [1 + °( a )] £k+i/2- This turns out to be inconsistent because then the con- 



dition (3.11) demands n to be half-integer which it can not be. If one takes the more general 
ansatz £k+i/2 = al £t+i/2 tne condition of uniform scaling of W( | p) yields the allowed sequence 



{n, 1} = {m + 1, 1/2}, {m, 3/2}, {m — 1,5/2}, . . .. The scaling of the lowest moments in eq. ( p.l2j ) 
remains unchanged, the uniform scaling however already starts with 5 n scaling like a 1 and thus 
simply reduces the number of double-scaling relevant terms. 
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the superloop equations (|2.11|) and (|2.13|) after a partial fraction decomposition. To 



optimize the procedure to only produce terms which are relevant in the double- 
scaling limit we have to analyze the operators appearing on the right hand sides of 
the superloop equations, i.e. the superloop insertion operators 5/5V(p) and 8/8^(p) 
as well as — Vo(p)). 

From the point of view of the a — > limit the effect of a given operator in 
8/8V(p) acting on an expression which scales with a to some power is to lower this 
power by a certain amount. Carefully examining each term in 8/8V(p) shows that 
a is maximally lowered by a power of (m + 3/2). All operators which do not lower 
a by this amount are subdominant in the scaling limit and may be neglected. The 
outcome of this analysis for 5/5V(p) is 

8 ~ 8M k 8 8x 8 



$V(p) x fi 6v (p)x 5M * SV(p) 8x 



where | 



+ ^- Al ) t + VA- J- (3.14) 
8V(p) x 8(E 1 - Ai) + 8V(p) x 8~ k 1 1 



m (k + b^ +1) (p) + (k + h^^Hp) (3.15) 



8V(p) x y 2^ x y 2 ' M 1 

8x 1 , , , 



TT^Kp) (3-16) 



8V(p) M 1 

" (3-17) 



^( S i ~ A = If* tO) 



8V(p) x 2 M x 



k 



(k-h^^ip) (3-18) 



and 



<t>M(p) = (p-x)- n - 1 / 2 d- 1 / 2 , (3.19) 



with d c = x c — y c . Repeating this analysis for the fermionic superloop insertion 
operator 8/8^(p) shows that here a is maximally lowered by a power of (m + 1) and 
the relevant contributions are 

* 6 f, 8E h 8 

8V(p) x SV(p) x 5 (Ei — Ai) ~ 5H* l ' ; 

with 



! Here the subscript x indicates that the critical behaviour is associated with the endpoint x 
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ggi - Ai 



-d c $\p) 
-d c ^t l \p)- 



(3.21) 
(3.22) 



Finally we state the double-scaling version of the operator (\l/ — t>o(p)) acting on 
the function <t> x n \p) 



($ - «b(p)) ^(p) 



E 

fe=i 



J n+2-fc 



4 (p - x)* 



+ 



AO 



2 4 



(P 



,n+l 



(3.23) 



Here the operator (^ — v Q (p)) is seen to increase the power of a of the expression it 
acts on by a factor of m. 

We are now in a position to calculate the double-scaling limit of the right hand 
sides of the loop equations (|2.11|) and (|2.13|) for u g (p) and v g {p) provided we know 

-i{p) and Viip), . . . , Vg-iip). As all the 



u. 



the double scaled versions of ui(p), 
y dependence has disappeared we do not have to perform a decomposition of the 
result. Moreover no and A*, dependent terms will contribute if we do not start 
out with any and we do not. 

The starting point of the iterative procedure are of course the genus correlators 
B. Keeping only the doublescaling relevant parts one has 



d_ 5v jg) 
dq S^(p) x 



1 



\p=q 



Suo(p) = 1 

6V{p) x ~ 8 ip-xf 



and 



Svp(p) 
5V(p) x 



(Hi - Ai) 



4d c Mi J ip 



+ 



x 



—2 



44M X J ip 



x 



(3.24) 



(3.25) 



The higher genus correlators u g (p) and v g (p) are expressed as linear combinations 
of the basis function X (p) an d take the general form B 



39-1 

E Af x {k) ip) 



3.y 



and v g ip) = £ X (fc) (p) + M (0) (p), (3-26) 



fc=i 



fc=i 



where k 9 is the zero mode coefficient not determined by the first two superloop 
equations Q2.11|) and ( |2.13| ). Note that the A g k ^ coefficients should up to a factor of 
two be identical to those of the hermitian matrix model obtained in ref . Q . We have 
calculated the A g k > and coefficients in the double-scaling limit for g = 1,2 and 
3 with the aid of Maple. The results of Ambj0rn et al. 
could be reproduced. 



for the A g k > coefficients 
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Before we state our explicit results let us turn to the general scaling behaviour 
of the one-superloop correlators 



W g (\p)~ a Q- 23 H m + 1 M- 1 , and ^(p|)~a (1 - 29)(m+1/2) - 1/2 , (3.27) 
which one proves by induction. The structure of the coefficients is seen to be 

B« = £ <«i, • • • , P I Mai M^ a 9 ' Ep , (3-28) 

a j >/3 1 c 

where the brackets denote rational numbers and where we write Hi for (Hi — A x ). 
One shows that the a, atj, j3 and s obey the conditions 

s 

a = 2g + s-l, and Y,( a i ~ X ) = 3 9 + 1 ~ f 3 ~ k ( 3 - 29 ) 

i=i 

with «j G [2, 3g>] and /? G [1, 3g]. For the zero mode coefficient n g the general struc- 
ture is given by the same expansion as eq. (|3.28 ) with k = 0. The conditions on a, 
aij, (3 and s then read 

s 

a = 2g + s, and £(a,- ~ l ) = 3 9 + 1 _ P ( 3 - 30 ) 

where aj G [2, 3g] and (3 G [1, 3(7 + 1]. Similar relations hold for the Ay^ [Q. 
The explicit results for the B)f> coefficients for g = 1 and g = 2 now read 



(i) 



(2) 



B? 



B 



(2) 



M 2 H 2 

+ „ „ , 9 , ) 



8M l( f c 8 Mi 2 d. 



Mo 



-Ai) 



16 Afrd c 
203 M 2 E 5 



+ 



8 M x 4 



145 M, 



Bp> = - 
105 H 6 



5 (Hi — A/ 



16 M 1 d c 
63 M 2 3 H 3 



128 4 Mi 128Mi & 4 128 4 Mi' 



+ 



32 Mi 6 4 2 



105 M 4 H 3 145 M3H4 

+ , 9 , , A + 



128 4 Mi 
87M 3 M 2 H 3 

2~ + 



+ 



128 4 z Mi 4 128 4 z Mi 
3 M 2 2 H 2 



105M 5 H 2 77 M 4 M 2 H 2 
' r ~ 32 Mi 5 d, 2 



75 M,M/ 2 H, 63 M 2 4 H 2 



63 M 9 2 H 



c 

2 ^4 



16 Mi 6 4 2 



32 Mi°4 
21M 2 3 E 2 



64 Mi b 4 2 64 4 2 Mi 



32M/4 2 



32 Mi 5 4 2 ' 



21M 2 H 4 77 M 3 M 2 H 2 
+ 128 Mi 5 4 2 



+ 



75 M>,.U 



2/ = 



Ai) 105Af B (Si-Ai 



+ 



32 Mi 6 4 2 
35 H 5 



35 M 4 H 2 
128 4 2 M/ 
63M 2 4 (H! 



Ax; 



256 d c 2 Mi 



128 M/4 



145M 3 z (Hi-A!) 21M 2 ^~ 3 



64M/4 2 
35M 3 H 3 



256 Mi°4 



64 Mi°4 128 4 Mi 



4 • 
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77 .U;.U,(E, A; 

5 J 2 



64Mi°d, 



c 

B ( 3 ) = 599M 3 M 2 (Si - Ai) lOSMa^-^Ai) 5M 3 ~ 2 

2 128Mi 5 4 2 32Mi 6 4 2 16 4 2 Mi' 
385M 4 (Si-Ai) 21M 2 2 S 2 7M 2 2 3 

H , 9 „ , /I "I" „ , ^ , 9 "I" 



256 4 z Mi 4 ' 64Mi 5 4 2 128d c 2 Mi 4 ' 
357Af 2 2 (Sx - Ai) 875 M 3 (Si - Ai) 352 3 63M 2 S 2 



(4) 

-D2 — „. rrm H 



2 64Mi 5 4 2 "'" 256rf c 2 Mx 4 128 M^ 2 128d c 2 Afi 4 ' 

(5) 105S 2 1617j^gi-Ai) (8) _ 1155 (Sj -AQ 

2 128M x 3 d c 2 256d c 2 Mi 4 ' 2 256 M^d c 2 1 ' J 

Note that the terms listed above are only potentially relevant, depending on which 
multi-critical model one wishes to consider. For an m'th multi-critical model all 
terms containing Mfc, k > m, or Hj, Z > m + 1, vanish in the double-scaling limit. 
We remind the reader that we assumed to have a non-symmetric potential and that 
the critical behaviour was associated with the endpoint x. In the case where the 
critical behaviour is associated with the endpoint y all formulas in this section still 
hold provided d c is replaced by — d c , by by A^ and x by y. 

3.3 The Iteration for F g , K g and u g (p) 

Having computed u g (p) and v g (p) (up to the zero mode) we may now proceed to 
calculate the free energy F g and the zero mode coefficient K g . This is done by 
rewriting u g and v g as total derivatives in the superloop insertion operators 5/5V(p) 
and 5/5ty(p) respectively, yielding the bosonic and doubly fermionic parts of F g as 
well as Kg. The procedure to compute the bosonic part of the free energy F g works 
just as in the hermitian matrix model described in ref. [[|. We have Fg OS = 2 Fg erm . 
From eq. ( |3.27| ) we see that F g scales as 

F g = F g bos + F g ferm ~ (2-«»)M/2) ) (3.32) 

just as the hermitian matrix model at its m'th multicritical point. 

To obtain the doubly femionic part F g m simply rewrite the basis functions 
^w(p) of v g (p) appearing in (|3.26| ) in terms of the functions <f>^(p) which by eqs. 



( p.21| ) and ( |3.22| ) are nothing but total derivatives in 8/5^f(p). Doing this for v g (p) 
allows one to directly deduce the form of _Fj erm and K g . 
The explicit results for genus one are 

pferm a \ / 5 5 4 3 M 2 S 3 5 M 3 S 2 3 M 2 2 ~ 2 



lQd^M^ ' 8d c a Mi a ' IQd^M^ %d c A Mf> 8 4 2 Mx 2 

(3.33) 



Actually the presentation here is somewhat misleading. In order to compute the u g (p) and 
v g (p) by iteration in the way outlined in the previous subsection one needs to know the full 
v\(p), . . . ,v g -i (p) including the zero mode coefficients k±, . . . , K g -i- In practice one hence computes 
all quantities at genus g, i.e. u g , v g , n g and F g , before proceeding to genus (g + 1). 



13 



and 

5 5 4 _ 3M 2 S 3 _ 5M 3 S 2 3M 2 2 S 2 
Kl 16M x 2 d c 8Mx 3 rf c 16M! 3 rf c ^' j 

This expression may of course be alternatively obtained by taking the double 
scaling limit of the genus one result of ref. computed away from the scaling limit. 

Before presenting explicit results for genus two, let us describe the general struc- 
ture of F s ferm for g > 1 

77>ferm \ ~> / n a I _,\ ^Pi ^/3 2 M ai . . . M aa /o or\ 

F g = 2^ ("l, ■ Pi, P 2 | Q)g a g+l , (3-35) 

where we write Si for (Si — Ai) and where the brackets denote rational numbers as 
before. The atj, Pi, a and s are subject to the constraints 

s 

a = 2g + s, and ^(aij - 1) = 3g + 2 - fa - fa, (3.36) 

i=i 

where oti G [2, 3g] and fej G [l,3j+ 1]. Similar relations hold for F^ os [|J. 
The result for genus two reads 



^fcrm 



(Hi - a x ; 



1015 E 5 M 3 375 H 4 M 2 M 3 
128 Mi 5 4 3 16 Mi 6 4 3 



315: 



M 2 C 



16 Mi 7 d c 3 
1155 H 7 
" 256 Mi 4 4 E 



385 H 4 M 4 1323 

+ \A + 



l 2 M 2 5 



64 Mi 5 4 3 



64 Mi 9 4 3 



693 E G M 2 
64 Mi 5 4 3 



63 E 2 M 2 M 5 525H 5 M 2 2 1155 H 3 M 2 M 4 



315 H a Ms 



4Mi 6 4 3 
2175 H 3 i 1 ^ 2 



32 Mi 6 4 



64MiX 



+ 



675 S 3 M 2 2 M 3 1785 H 2 M 3 M 4 



64Mi 5 d c 3 256Mi 6 rf c 3 ' 16Mi 7 rf c 3 128Mi 6 4 3 
495 S 2 M 2 3 M 



8Mi 8 4 3 
1155 S 2 M 6 



/:, 2205 S 2 M 2 2 M 4 



1323 H 3 M 2 4 8175 E 2 M 3 2 M 2 



256 Mi 5 4 



64Mi Y 4 3 64Mi 8 4 3 ' 256Mi 7 4 3 

21 S 2 S 5 M 2 105 H 2 S 6 145 E 2 E A M 3 



16Mi & 4 



105~ 2 ~ 4 M 2 2 21 E 3 S 4 M 2 

"1" _ . - - R - 1 



128 Mi 4 4 3 
35 s 3 



128 Mi 5 4 3 
63 E 2 E 3 M 2 3 



64Mi 6 4 3 
35 E 2 E 3 M 4 



32 M 1 b d c 



+ 



64Mi 5 4 3 128Mi 4 4 3 32Mi 7 4 
195 E 2 E 3 M 2 M 3 



64 Mi d c 



(3.37) 



Of course so far we have determined only the coefficients F g of the genus expan- 
sion of the free energy (cf. eq. (|2.7| )). For an m'th multicritical model the rele- 
vant expansion parameter in the double-scaling limit is the string coupling constant 
A^~ 2 . If we introduce the bosonic and fermionic scaling moments \x\~ and 



a 



-2m-l 



T k by (cf. eqs. flp) and Q) 
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M k = a m ~ k ^ ke[l,m], ( 3 - 38 ) 
(Si - Ax) = a m+1/2 n, E t = a rn - l+3/2 r h I G [2, m + 1], 

we get by replacing Mk by /i^ and setting Mk equal to zero for k > m as well as 
replacing by t\ and setting S; to zero for / > m + 1 in the formulas above exactly 
the coefficients of the expansion in the string coupling constant. Needless to say 
that these results apply for non-symmetric potentials where the critical behaviour 
is associated with the endpoint y as well by performing the usual replacements. 
For the zero mode coefficient at genus two we find 



495 H 2 M 3 M 2 3 63 H 2 M 5 M 2 315S 3 M 5 1323 S 3 M 2 4 
K2 ~ 8 d 2 M^ + Ad 2 M^ ~ 64 4 2 Mi 5 + QAd^M^ 

675S 3 M 3 M 2 2 2175 S3M3 2 8175 E 2 M 3 2 M 2 385S 4 M 4 



16c/ c 2 M! 7 256 d 2 M^ 256 4 2 Mx 7 64d c 2 Mi f 

1 / _ 



2205 H 2 M 4 M 2 2 375H4M3M2 1155 S 2 M 6 1785S 2 M 4 M 



~t~ .. „ ■ 9 , ft „_„,9,.-^~t~ 



64 4 2 Mx 7 16 4 2 M! 6 256 d c 2 Mi 5 128d c 2 Mi 

693 5 6 M 2 315 5 4 M 2 3 525H 5 M 2 2 1015 H 5 M 3 



64 d/MS \§d 2 M x l 32d c 2 Mx b 128 
1323 S 2 M 2 5 1155S 3 M 4 M 2 1155 5 7 



64 4 z M! y + 64 4 2 Mx 6 + 256 4 2 M! 4 ' ^ 3 ' 39 ^ 

We obtained these results with the aid of a Maple program which performs the 
iteration up to arbitrary genus. | In practice the expressions become quite lengthy, 
e.g. _F| erm consists of 114 terms. 

For the knowledge of the full W g ( \ p) it remains to compute u g (p). This is of 
course done by applying 5/5V(p) to Fj erm . The general structure of u g (p) turns out 
to be @ 

3<?+l 

«»= E^'^W, (3-40) 
fc=i 



where 



underlying the conditions 



Af = £ ...,«., A, ft 1 a),,» Mai fffi (3-41) 



a = s + 3#, and - 1) = 4 + 3g - k - ft - ft, (3.42) 

i=i 

with OLj G [2, 3(7 + 1] and ft G [1, 3g + 2]. Due to space let us only state the genus 
one results 



5 The reader interested in higher genus results should feel free to contact the author. 
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= - 

1 32d c 2 M! 2 16d c 2 Afi 3 32 4 2 M! 2 

lSMafSi-AOSg , 15M 2 (Si - Ax) H 4 9 M 2 2 (H 1 - A x ) S 3 



5S 2 S 4 3M 2 H 2 H 3 35(Si-Ai)S E 



~ {2) 15(Si-Ai)»4 3S 2 H 3 3M 2 (Si - Ai) » 3 

1 16 4 2 Mi 2 8 4 2 Mi 2 44 2 Mi 3 

5^,(5,-^)5, _ 3Af,'(5,-A,)5, 

~ (3) ^(Si-AQSa 25M 2 (S 1 -A 1 )S 2 ^ (4) = SS^-AQS, 

1 16 4 2 Mi 2 32 4 2 Mi 3 ' 1 ~ 32 4 2 Mi 2 

This concludes our analysis of the double-scaling limit for generic potentials. 

The case of symmetric bosonic and generic fermionic potentials was considered 
in ref. [fjj where a doubling of degrees of freedom was observed for genus zero. With 
the methods presented in this paper one can see that this holds for higher genera as 
well, i.e. the free energy here takes the form 

F symm = 2 pbos + ^r m ,x + ^fernx.y ^ (3 44) 

where F^ os denotes the bosonic part of the free energy for generic potentials, i^ erm > x 
and Fg m ' y denote the doubly fermionic parts in the generic case where the critical 
behaviour is associated with the endpoints x and y = —x respectively (cf. eq. ( |3.35| )). 
If one chooses to take the fermionic potential to be symmetric as well the doubly 
fermionic part of the free energy will vanish, which may be directly deduced from 
the results of ref. J?]. 



4 Conclusions 

We have studied the double scaled supereigenvalue model in the moment descrip- 
tion. The m'th multi-critical point was identified and the scaling properties of the 
moments and basis functions derived. The iterative scheme for the calculation of 
higher-genus contributions could be optimized to produce only terms relevant in 
the double-scaling limit. The general form of the free energy and the one-superloop 
correlators at genus g were found. We presented explicit results up to genus two. 

We believe that this paper shows once more the effectiveness of the iterative 
scheme and the moment description of Ambj0rn et al. The analogy of structures 
in the hermitian one matrix model and in the supereigenvalue model continues to 
hold. It is interesting to note that the analogous computations for the hermitian 
matrix model [HI were the basis of a transition from the double scaled one-matrix 



to the Kontsevich model |T2[]. One might speculate that our results are connected 



to properties of the moduli space of super-Riemann surfaces. 
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There are more interesting unanswered questions to be adressed, such as the 
supersymmetric generalization of two and multi-matrix models, where one would 
hope to reach new minimal superconformal models, or the generalization of matrix 
models in external fields, perhaps a step towards a matrix based formulation of the 
supereigenvalue model. 
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